ON THE BOUNDEDNESS OF AN ITERATION INVOLVING 
POINTS ON THE HYPERSPHERE 

THOMAS BINDER AND THOMAS MARTINETZ 

Abstract. For a finite set of points X on the unit hypersphere in K d we 
consider the iteration iij+i — Ui + x%, where Xi 1S the point of X farthest from 
itj. Restricting to the case where the origin is contained in the convex hull of 
X we study the maximal length of U{. We give sharp upper bounds for the 
length of Ui independently of X. Precisely, this upper bound is infinity for 
d > 3 and y/2 for d = 2. 



1. Introduction and overview 

Throughout this paper we will assume that d > 2. By M. d we denote (i-dimensional 
Euclidean space, equipped with the standard scalar product (-, •) and induced 
norm || ■ ||. Moreover S l (r) denotes the /-dimensional sphere of radius r, and 
S l := 5^(1). These spheres are always considered as embedded in M. d . Let 
X = {x\, . . . ,x n } C S^ -1 C M. d be a finite set on the unit hypersphere. With- 
out mentioning this each time, we assume that the linear space spanned by the 
elements of X equals M. d , i.e. d cannot be reduced. Consider the iteration 

u := 0, u i+l := Ui + Xi, 

where % G No and Xi is the element of X which is farthest away from Ui (which 
happens to be argmin a;gX (x, ui)). In case there are several elements of X at 
maximal distance, just choose any of them. Due to this ambiguity there are 
many iterations («j)°i for a particular set X. By U(X) we denote the set of 
vectors occurring in any of these iterations. Let 

u*(X) := sup { | u G U(X)} 

be the greatest length reached during any of these iterations. The question which 
values u*(X) can take is simple and intriguing; it was brought up in connection 
with the rate of convergence of an iterative approach of computing the smallest 
enclosing ball of a point set, as described in the following. 

Let Y C M. d be a finite set of points. Then the smallest enclosing ball SEB(y) of 
Y exists and is unique |Wel91] . We assume that Y has at least two elements. By 
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cGl d and R G M + we denote center and radius of SEB(F), respectively. Badoiu 
and Clarkson |BC03j introduced the following approximation of c: 

c := 0, Cj+i := q + — — (& - Cj), (1) 

z + 1 

where i e N and is the element of Y farthest away from Cj. This approximation 
(cj)i^o ^ s related to the iteration (ui)°Z by Rui = z(q — c) which implies = 
Ui + The set X connected to (ui)°Z is given by 

X:={hy-c) \yEY). (2) 



R 

Unlike X the set X can contain also points in the interior of the unit hypersphere. 
Martinetz, Madany and Mota [MMM06J show that after a finite number of steps 
all £j will lie on the boundary of SEB(Y), i.e. ^ gF for all i > io, where Y C Y 
consists of all points on the surface of SEB(Y). This clarifies the correspondence. 

While the approximation is extremely easy to use, the question of convergence 
needs to be answered. In [BC03] it is shown that for i G N 

lie — cv II 1 

-IT 1 £ Vi (3) 

[MMM06J aims at proving faster convergence than ([3]). In particular: 

Theorem 1 ( |MMM06| . Theorem 2). Let Y CR d be a finite set with at least two 
elements, and let X be given by (J2J). Consider the approximation (TjQ) o/SEB(Y). 
Then for all i G N 

II c — Ci\\ u*(X) 



R ~ i ' 

where the definition of u* has been extended to sets X with points on or in the 
interior of the unit hypersphere in a straightforward manner. 

In view of Theorem [H a finite value of u* or even a uniform upper bound inde- 
pendent of X is desirable. Before stating our results on the latter, we need some 
preparations. 

The connection between (cj)^ and (wi)^Q is further illustrated by 

Proposition 2. For a finite set X C S' 1 ' 1 C IR d the following statements are 
equivalent. 

(i) SEB(X) = S d -\ 

(ii) The origin G M. d is contained in conv(X), 

(iii) 5(X) > 0, where 

S(X) := — max min(x, u). 

\\u\\=ixex" 
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Proof. (i)<^=^(ii) is due to R. Seidel (cf. Lemma 1 in [FGK03J). (ii)^=^(iii) 
follows from the fact that a point p G W d lies in the convex hull of X if and only 
if min^gx (x — p,u) < for all unit vectors u. □ 

X is called O-balanced if ^ conv(A). For 1 < b < d — 1 the set X is called 
6-balanced, if is a point on the boundary of conv(A) and is contained in a 
^-dimensional face, but not in a (6 — l)-dimensional face of conv(A). If is an 
inner point of conv(A), then X is called (/-balanced or balanced. Having the same 
balance property is an equivalence relation on all sets X under consideration. 

Note that 5(X) is strictly positive if and only if X is (i-balanced, and Proposition [2] 
characterizes all sets X that are not O-balanced. 

Theorem 3. Let X be a finite set of unit vectors in M. d . 

(i) If X is O-balanced, then u*(X) = oo. 

(ii) If X is b-balanced for < b < d, then u*(X) < oo. 

Proof. Again, (ii) is shown in [MMM06J; it remains to prove (i). Since conv(X) 
is compact, there is a point T £ conv(X) which is closest to the origin. Let 
e := \OT\. Clearly \\xj\\ > e for all j G N , therefore = || 5^y=o Xi 1 1 — ^ e * s 
an unbounded sequence for i e N . □ 

For < b < d we define 

u* d * b := sup { u*(X) | X C S*- 1 C R d finite and fo-balanced}. 
Our goal is to compute u* d * b for all possible d and b. 

Theorem 4. For d = 2 we have u^q = oo, while = = 

Clearly, for d = 2, X = {xi,x 2 }, X\ = (0,1), x 2 = (1,0) the iteration u = 0, 
Ui — Xi, U2 — %i + x 2 is valid and ||m 2 || = y/2. This manifest example represents 
one inequality of the proof of Theorem HI the missing inequality is shown in 
Section [2j 

Theorem 5. For d > 3 we have u* d * b = oo for all < b < d. 

Proof. For any dimension d we have u* d * = oo from Theorem [3] (i). For 1 < b < 
d — 2 the assertion follows from the example discussed in Proposition [13] below. 
For b = d and b = d — 1 use Proposition [T^] (ii) and (iii), respectively. □ 

Although the balance property of A is a suggesting geometric property, it does 
not seem to give a finer prediction for u*(X) than o~{X). In the balanced case, 
< S(X) determines a finite upper bound for u*(X) as shown in [MMM06J, 
namely 

IWI £ Wx) + 1 ' 8 e No ' 
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With respect to the faster convergence we have an immediate result for d = 2: 

Corollary 6. Let Y C M? be a finite set with at least two elements. Assume that 
all elements of Y lie on the boundary o/SEB(F). Then \\c — q|| < for all 
i G N. 

2. Proof for d = 2 

Let ei, e2 denote the canonical orthonormal basis of R 2 . Each Xj 6 X, 1 < j < n 
can be written as 

xj = cos(0j) ei + sin(0j) e 2 = [1; 

where [f; 0] indicates a point in standard polar coordinates on R 2 . Similarly, for 
j G N we write 

Xj = cos(%) ei + sin(^) e 2 = [1; tpj], 

Uj = Aj(cos(aj) e\ + sin(aj) = [Xj] a.j\. 

All argument angles are real numbers taken modulo 2n. The freedom in rotation is 
fixed as follows. Assume that numbered counterclockwise, starting 

at 0i = 27r — 0, ending at <ft n = n + <fi, such that there is a gap with angle size 
7r — 2(p between the two neighboring elements x\, x n of X is symmetric about 
the e2-axis. We call this a parametrization of X with base gap of size tx — 2(f), 
where (f> G [0, |). The choice of (j) indicates that we restrict to the balanced cases. 
Define (f> := f - (j). For C M 2 and fc = 1, . . . , n let 7fc(W) denote the set 
obtained by translation of W by x^. The set T is defined by 

T:= {[f;0] GM 2 | fe(l,V^] and G (f ~ 0, f + ^)}- 
Moreover, we define three subsets of M 2 by 

R := {[?] <f>] I f > and </> G (tt - 0, 2tt + 0)}, 
Q := {{a, b) I |a| tan0 < b < \a\ tan0 + X m i n }, 
P:= {tiGM 2 < 1}\(RUQ). 

Here A m j„ := 2< ^^ is the length of the intersection of Q with the e2-axis. Figured] 
gives an illustration of this situation; |FIG] gives an animated version where 
varies in time. 

Lemma 7. Let X be a finite subset of S 1 C R 2 , parametrized as above. Suppose 
that G [0, |), ie. t/ie size of the base gap is greater than |7T. Define the set V 
by 

V:=PU T n (P + ) U Ti(P-) UQUi?, 
where P + , P~ denote the elements of P with non-negative and non-positive ex- 
coordinate, respectively. Then Uj G V for all j G Nq. 
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Figure 1. An arbitrary set X C S 1 C IR 2 given in base gap 
parametrization. Only x\ and displayed, the remaining el- 

ements of X are above x\ and x n . Recall that <f> + <fi = |. R is 
the open set bounded from above by the lower dashed lines. Q is 
the closed set between the dashed lines. The set P is given by the 
central hatched area. For small values of (ft, 1\{P~) \{Q U R) and 
T n (P + ) \(Q U R) are nonempty. 
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Proof. Clearly ito G V. By induction, assume that Uj G for some j G N. The 
proof is complete if all of the following claims are shown to be true. 

(a) If Uj G Q, then -u J+1 e Q U R. 

(b) If Uj G P, then G T„(P + ) U Ti(P~). 

(c) If Uj G P, then GPUQUi?. 

(d) If uj G T n {P + ), then G P U Q U P. 

(e) If Uj G Ti(P'), then G P U Q U R. 

If Uj £ P U Q, then x x or x n is chosen in the next step of the iteration, i.e. 
Xj G {xi,x n }. Therefore, (b) is trivial. Also (a) is true since 7i(Q) and T n (Q) 
have no parts above Q. If (d) is true then (e) holds by symmetry. Hence it 
suffices to show (c) and (d). 
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Claim (c). Suppose that Uj G R is arbitrarily fixed. If aj G <p, 2tt — </>), then 
from Figure [T] it is clear that translation of the part of R with such argument a,j 
by an arbitrary unit vector stays inside P U Q U R. 

Otherwise, etj G [—0, 0) or dj G (-7T — 0, ir + 0] , where the second part follows from 
the first by symmetry Restricting to a := ctj G [—0, 4>) and setting A := \j > 0, 
•0 := ipj G [7T + 2a — <fi, 7r + <fi] we can write 

Uj+i = (A cos a + cos^)ei + (A sin a + sin-i/*^. 

The range of ip follows since the center of the interval of possible values for ifi is 
a + 7r, it extends by % + (/)— (a + n) = <p — a to both sides. We continue to work 
on two cases. 

(c.i) The ei-coordinate of Uj + \ is non-negative. In this case sin(-0 — <p) < ^ 
and Asin(</> — a) > 0. Since equality does not hold simultaneously, 

< A sin(0 — a) + sin(0 — ip) H — — . 

Expanding and rearranging the trigonometric terms, substituting A m j n = 
(which denotes the length of the intersection of Q with the e2-axis) 
and dividing by cos <fi > we get 

(A sin a + sin tjj) — X min < tan (A cos a + cos ip) . 

This shows that Uj + i falls below the line bounding Q from above. Hence 
u j+1 G Q U R. 

(c.ii) The ei-coordinate of Uj + i is negative, i.e. A < — If we knew the 

inequality 

COS-0 , 

->2cos^-a, 4 

cos a; 

then A < — 2 cos(-0 — a) would follow using the inequality for A. We would 
arrive at 

||M i+ i|| 2 = 1 + A 2 + 2Acos(> - a) < 1, 

which would show that Uj + \ G P U Q U R. Hence we are left with (J4]). 
First consider the case a > 0. Then 2 cos(-0 — a) < —a/3 and 

cos -0 12 
cos a cos a 3 

hence (jl]) is true for this case. Now restrict to the case when a < 0. Then 
2 cos('0 — a) < — 1 and 

cos-0 ^ cos(-7r + 2a — (f)) 
cos a — cos a 

hence (HI) is true. 
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Claim (d). From the assumption there is some v = [A; 5] G P + with 2s hi(f- < ft) — 
A < 1 and 5 G [f - 0, f ] such that 

«j = = (A cos 5 — cos0)ei + (A sin 5 — sin0)e2- 

We are done if we show that Xi is chosen for the next step of the iteration, i.e. 
Xj = %i- In this case 

Uj +1 = A cos 5ei + (A sin 5 — 2 sin 0)e 2 . 

uj + i has a smaller e2-coordinate than the original point v G P + , hence Uj+i G 
iiUQU P + . We are left with the mentioned claim and show that the argument 
angle ctj of Uj satisfies a>j < tt — 0. From 

, . rx v3 sin(0 + 5) \/3 . „. 

Asm + 5 > ^- . ; > ^- > sm2</) 
2 sm(o — 0) 2 

we get 

(A cos 5 — cos 4>) sin > — cos 0(A sin 6 — sin 0) . 

Since A sin 5 — sin <p > and sin > division by these terms does not change 
the type of inequality. We obtain 

A cos 5 — cos (j) . / , \ 

COt Otj = > — COt = COt(7T — 0), 

A sin o — sin 

which proves the desired fact. □ 
Lemma 8. In the situation of Lemma\^ we have V fl T = 0. 

Proof. By construction (PUQUR)nT = 0. By symmetry it is therefore enough 
to show that T n {P + )nT = 0. As before, let u = [A; 5] G P + , where S G [f - 0, f ] 

and ^#-7, < A < 1. Then 

2sin(d— <p) — — 

T n u = (A cos S — cos0)ei + (A sin S — sin0)e 2 . 

Starting with 

^3 

Acos(<5 — 0) < cos(£ — 0) < < cos(0 — 0), 
expanding and dividing by A sin 5 — sin > and by cos > we get 

A COS 5 — COS - / 7T 

cot arg T n u = < — tan = cot — h 

A sin o — sin V 2 

which shows that the argument angle of T n u is greater or equal than | + 0. 
Therefore G" T, which proves the assertion. □ 
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Proof of Theorem^ Again, the set A 2 ^ from Example [TOl below shows that u 2 \ > 
\/2. Moving t\ slightly away from e 2 turns A 2j i into a balanced set and shows 
that also u* 2 * 2 > v2. Hence it suffices to prove u 2 \ , u 2 * 2 < Contrarily, we 

assume that there exists an iteration such that A^ > y/2 for some fixed i G N. 
Without loss of generality we may assume that i is the smallest such index, in 
particular Aj_i < v2. 

The angle jj G [0, 7r] between Uj and Xj is defined for all j G N since without loss 
of generality we may assume Uj 7^ 0. Now observe that 

! + 0=i(27r-(7r-20))< Ti <7r 

for all j G N. A simple computation yields 

A 2 = 1 + 2Aj_i cos 7j-_! + X)_ v (5) 



Hence 
and 



2Ai_i cos 7^1 = A 2 - A 2 _ x — 1>2 — 2 — 1 = — 1, 



11 /7T 

-- < "T^T - < cos 7i-i < COf %2 + ^' = ~ slIli 



since from §5§ we also have 1 < Aj_i. Therefore 

7T 2 7T 

— h © < 7i_i < -7r and < < — . 

In other words there is a gap greater than ~7r between two neighboring elements of 
X. In a second step of the proof we will explore possible ranges of a^i. Clearly, 
the angle between Ui-\ and X ^ - X iyi IS less or equal than ~7r. Therefore exactly one 
of the following cases holds. 

Case 1. G (| — 0, I + 0), where := | — 0. Hence G T but also 

G from Lemma [71 This contradicts Lemma [H 

Case 2. oti-i G (|7r — 0, ~7r + 0), where := | + 0. We can restrict the range 
of aj_i further by adding the above condition not only for X\ and x n , but for all 
elements of X. Doing so we get that 

|7r > — <pj, if 7r > — <pj, and 
< cti-± - (f)j, if 7r < - (f)j. 

Let = l,...,n- 1 be the greatest index satisfying n < — 0^. Since k is 
maximal we have tv > — <fik+i- We get <fik+i — <fik > f 71 ") which shows that 
there must be a second gap which is greater than |7r. After a rotation of the 
coordinate system and renumbering the elements of X we may apply Lemma [S] 
again and obtain a contradiction. 

The indirect assumption must have been wrong in Cases 1 and 2, hence both 



u 



2,1> "2.2 



< V2. □ 
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3. Examples 

This section provides examples illustrating that the situation is more complicated 
in dimension d > 3. All examples are unique up to rotation of M. d . 

Example 9. For I > 1 we describe the operation of choosing / + 1 equidistant 
points xq, . . . ,xi E S 1 ^ 1 C W. Equidistant means that the value s of the scalar 
product does not depend on the chosen pair of points. Since all vectors have unit 
length, the constant scalar product equals cos a for some a G [0,7r]. By recursion 
on I suppose xi, . . . , x\ have been found in the next lower dimension I — 1, with 
scalar product s. Set 

Xq — (0, 0, . . . , 0, 1), X\ = {x\ cos a, sin a), xi = (xi cos a, sin a). 

We demand 

sin a = (xo, x±) = s = (x,, Xj) = sin 2 a + (xi, Xj) cos 2 a, 

which leads to s = s 2 + (1 — s 2 )s. Solving this equation gives s = j^. It is easy 
to see that the recursion produces the values 

-1 -- -- -- 

' 2' 3' 4'"" 

for s. Hence, when denoting the scalar product of dimension I by s/, we get 
si = s = — j. Knowing s it is also clear that xo + . . .+Xd = since x%+. . .+Xd = 0. 
In low dimensions, equidistant points are just two points on the real line (/ = 1), 
a regular triangle in a circle (I = 2), or a tetrahedron in a 2-sphere (/ = 3). 

Clearly, the set X of d + 1 equidistant points is balanced in S^ 1 C M. d . The 
problem of finding u*(X) in this case was approached by a computer experiment 
only. We checked d = 2, ... ,12 and found that u*(X) = ^p, where a is the 
integer sequence 

0, 1, 2, 4, 6, 9, 12, 16, 20, 25, 30, 36, 42, . . . 

starting at index d = 0. Obviously, itj may take only a certain finite number of 
values on the lattice 

d+l 

{^hxi | fcj eNo}, 

1=1 

all of which are close to the origin. For example, there are 3 possibilities for d — 1 
and 7 for d = 2. The sequence a has relations to other fields and problems p\TTj. 
Note also that a(d) < dyfd, or equivalently u*(X) < y/d. The latter inequality 
was an ad-hoc conjecture for a general set X, which turned out to be true only 
in dimension d = 2. 
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Example 10. For 1 < m < d consider the following set X = t7n consisting of 
n = d + m points. As before, let G M, d be the vector with all zero components 
except the ith which is 1. Then define 

Ad, m ■= {ei, e 2 , • • • , e d , —e\, -e 2 , • • • , -e m }. 

Proposition 11. Let X = A^ m be as in Example ITU 

(i) Ad, m is m-balanced, 

(ii) u*(A djTn ) > y/d-m+1. 

Proof, (i) is clear from the definition; the origin is contained in the m-dimensional 
face of conv(A rfjm ) spanned by ±ei, . . . , ±e m . For (ii) observe that there is an 
iteration such that iij = e m+ i + e m+ 2 + . . . + e m+ i for 1 < % < d — m. □ 

It is likely that equality holds in (ii), but we do not need this stronger assertion. 

Example 12. The following construction of X = B^ ^e, 0) depends on the di- 
mension d, some integer 1 < b < d — 2, some real numbers e > and < < ~, 
where the value of is uncritical. For c:= d — b, 2 < c < d — 1, we have the or- 
thogonal decomposition IR d = M. b © M c . The subspaces contain unit hyperspheres 
S^ 1 C and S^ 1 C| c . 

In S^ -1 choose c+1 points follows. Fix any direction v e S c 1 and 

consider the linear hyperplane V which is perpendicular to v. In S c ~ 2 = VClS^ 1 
choose c equidistant points described in Example [HJ Then let 

Xi := cos(e) Xi + sin(e) i> 

for % = 1, . . . , c. Note that equidistant in S c 2 (cose) := (V + 

sin(e)v) fl if? 0-1 . The remaining point xo is given by 

x := — cos(0) Xi + sin(0) v. 

In S" 6-1 choose 6+1 equidistant points x c+ i, . . . , Xj+i, which makes a total of 
n = d + 2 points in X. 

Proposition 13. For d > 3 and X = B^ ^e, 0) £/ie following statements are 
true. 

(i) X is b-balanced, 

(ii) /or any /arge M > there is an e > sttc/i iaat n*(X) > \[M. 

Proof, (i) is clear from the definition; the origin is contained in the o-dimensional 
face spanned by x c+ i, . . . , Xd+i- Note that x\ + . . . + x c — csin(e)t> and 

a := (xi, Xj) = (xi, Xj) cos 2 e + sin 2 e = 1 cos 2 e 
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since (xi,Xj) = — for all 1 < i,j < c. From now on we suppose that e is 
sufficiently small such that 

- -— < a < 0. (6) 

We also have 

-cos0 + sin (f> sine; i = 1 , 
—<7 cos (j) + sin sine; 1 < i < c. 



(x ,Xi) 



To prove (ii), we show that the iteration which starts with xq and adds points 
from {x\, . . . ,x c } as long as possible is feasible. More precisely, 

U = 0, Ui = Xq, U 2 = X + Xi, U C+1 = X + X 1 + h X c . 

In general for i = 0,1, . . . we can write 

Uic+i = x + (i - l)(xi + x 2 + ... + x c ), 

U ic+2 = Xq + (i - l)(xi + X 2 + ■ ■ ■ + X c ) + Xi, 

Uic+c = x + (i - l)(xi + x 2 + ■ . . + x c ) + (x 1 + x 2 + . . . + X c _i), 
M( i+ i) c+ i = x + i{xi + x 2 + . . . + x c ). (7) 

In what follows we fix < % < k and < j < c — 1 arbitrarily, and consider step 
s := {i + l)c + j + 1 of the iteration ([7]). In other words, we want to control the 
iteration up to and including step (k + l)c + m + 1, where < m < c — 1. 

(a) To be able to choose in step s we must have 

(u s ,x j+1 ) < 0. 

(b) Also, to make the choice of Xj + i work, the scalar product with all other 
vectors must be at least as big as the one from (a), or 

(u s ,xi +1 ) > (u s ,x j+1 ) 

for all < I < c - 1. 

(c) The point xq must not come into play, which is the case when 

(u s ,x ) > 0. 

(d) By construction we have 

(u s ,x r+1 ) = 

for c < r < d. 

Let us now analyze these conditions. There is nothing to show for (d). For (c) 
we compute 

, , _ J 1 + zcsin esin 0; j = 0, 

\u s ,xq) < j — cos (j) -\- ic sin e sin <p — (j — 1)^ cos (ft ~\~ j sin e sin <f>', 0<j<c— 1. 
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From this expression it is clear that (c) is always satisfied. Looking at (a) and 
(b) and observing that 1 + (c — l)cr = csin 2 e we compute 

, \ — j zcsin 2 e — cos0 + sin sine; j = 0, 

s ' J+ ' 1 i csin 2 e + jcr — a cos + sin sine; 0<j<c— 1 

and for I ^ j 

{icsin 2 e +(j — l)a + 1 — cos + sin sine; = I < j, 

icsin 2 e +(j — l)a + 1 — a cos + sin sine; < I < j, 

icsin 2 e +jcr — a cos + sin sine; I > j. 

From these expressions (b) is immediately clear; one just has to compare the 
varying terms and to use (Q. It remains to analyze Condition (a). For j = it 
can be expressed as 

cos — sin sin e , N 

i < 



• 2 ' V J 

csin e 

for j > note that we have a set of c — 1 inequalities, whose "sharpness" increases 
with j, cf. (jSJ). Therefore it suffices to take the last condition (J — c— 1) which 
reads 

cr(cos0 — (c — 1)) — sin sine 

% < — 2 . (9) 

csin e 

In the second and last part of the proof, the assertion is brought into play. Assume 
the length \fM is reached in step (k + l)c + m + 1, i.e. 

\\u {k+1)c+m+1 \\ 2 > M. (10) 
For arbitrary k and 1 < m < c - 1 we have 

H^O+ijc+m+iH 2 = 1 + (kc + 2m)kcsin 2 e + (l + (m — l)cr) (m — 2 cos 0) + 
2(kc + m) sin e sin 0, 
while for m = we get the simpler expression 

\\u(k + i) c+ i\\ =l + kc sin e + 2A;c sin e sin 0. (11) 

Assuming m = (to use the advantages of the simpler form) and inserting (fTTl) 
into (TlOl) we get an inequality which is quadratic in k: 

,, , 2 sin 1 - M 

fc + A; — — + -= 5- > 0. 

c sin e c 2 sin e 

Solving the inequality gives 

fe > V sin 2 — 1 + M — sin 
c sin e 

To finish the proof, we must put together (jSJ) and ffl2|) as well as (jHJ) and f|T2|) . 
For the first pairing, solve 



2 * cos — sin sin e 

sin — 1 + M — sin < . 

sin e 
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Isolating M yields 

1 



M < cos 2 0(1 + 



sin 2 e 



For small e, the right-hand side becomes arbitrarily large, which finishes this part 
of the proof. For the remaining pairing, one has to solve 



sin 2 (ft — 1 + M — sin < 



<t(cos0 — (c — 1)) — sin (p sine 



sm e 

Isolating M again gives 

a 2 (cosd> — (c — l)) 2 9 , 

M < — — V — + cos 2 </>, 

sin e 

which with small e again has an arbitrarily large right-hand side. □ 

Example 14. The following construction of a point set X = Cd(e, fi, 4>) depends 
on the dimension d > 3, on real numbers e > 0, fi > and < </> < |, where 
the value of 4> is uncritical. Pick any unit vector v G M d which determines a 
hyperplane V of M. d . In S" 1 ' -2 C V choose d equidistant points 
described in Example [9j Then define 

Xi := cos(e)xj — sin(e) v 

for z = 1 , . . . , d. The two remaining points are given by 

Xd+i = — cos(/i)xi + sin(//) v, 
xo = cos(0)xi + sin(0) v. 

Finally let X := {x , X!,...,x d , x d+1 }. 

Proposition 15. For d > 3 the following statements are true. 

(i) Cd{e, ix, (p) is d-balanced for e > 0, and (d — l)-balanced for e = 0, 

(ii) for any large M > there is an e > such that u*(Cd(e, 3e, |)) > \[M, 

(iii) for any large M > t/iere is a /i > swc/j £/ia£ M*(Cd(0, /i, |)) > \/M. 



Proof, (i) is immediately clear from the definition, in particular for e = the 
origin is contained in the (d — l)-dimensional face spanned by x±, . . . , x d . We are 
left with (ii) and (iii) which are shown simultaneously. Consider the following 
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finite piece of an iteration for C d (e,fi, (/)). Start with u = 0, and let 

Ul = X , 

U 2 = Xq + X d+1 , 

U-i = Xq + Xi + X d+1 , 

U2k-1 = X + (k - 1)(X 1 + X d+1 ), 

U2k = x + (k- l)(xi + X d+ i) + X d+ i, 
u 2 k+\ = x + k(x x + x d+ i). 

The following conditions (a)-(c) are sufficient for the iteration to work as above, 
up to step 2k + 1. 

(a) We must have (ui,xq) > for all 1 < I < 2k + 1, i.e. xq is never chosen 
between steps 2 and 2k + 1 of the iteration. 

(b) Additionally, also the scalar product with the other vector must be at least 
as big as the chosen one, meaning 

(u 2 i,xi) < (u 2i ,x d+1 ), (u 2i+1 ,x d+1 ) < (u 2i+l ,xi) 

for all 1 < % < k. 

(c) To be able to choose x d+ \ in step 2% and x\ in step 2% + 1 we must have 

(u 2h xi) < (u 2i ,x m ), (u 2i+1 ,x d+1 ) < (u 2i+1 ,x m ), 

for all 1 < % < k and 2 < m < d. 

In order to examine Condition (a) it is straightforward to compute 

, \_f 1 — cos(0 + e) + i(cos(0 + e) — cos(0 + /i)); I = 2i, 
{u l ,x )-^ 1 +i(cos(0 + e)-cos(<^ + /i)); l = 2i + l. 

Since /i > e for both (ii) and (iii), the terms on the right-hand side are always 
non-negative. Therefore (a) does not impose any additional condition. Similarly, 
for Condition (b) we compute 



{u u xi) 



cos(0 + e) — 1 + i(l — cos(yU — e)); I = 2i, 
cos(0 + e) + i(l — cos(yU — e)); l — 2i + l, 



I \ — / cos (/ i — e ) — cos(0 + /i) + i(l — cos(fi — e)) ; / = 2i, 

{Ul,Xd+l) ~ \ -cos(4> + fi) +i(l-cos(^-e)); Z = 2i + 1, 

which is equivalent to 

cos(0 + e) — 1 < cos(/i — e) — cos(0 + ^), 
— cos(0 + fi) < cos(0 + e). 
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Again, since both inequalities are always true, (b) does not introduce new con- 
ditions either. Finally, Condition (c) requires 

(ii 2 i, x m ) — (u 2 i, Xi) — — cos e(cos e — cos + z(cos/i — cos e)) > 0, 

a — 1 

(u 2 i+i,x m ) - (u 2i+ i,x d+1 ) = - - - - cos cos e + cos(0 + e) + cos(0 + ^) + 

d 

i- — -(cos/x — cos e) cos e > 0. 

We demand that if i satisfies the first inequality, then it shall also satisfy the 
second. This leads to the additional condition 

cos - cos e cos cos e - cos(0 + e) — cos(0 + y) 

cos \i — cos e ~~ ^-(cos fi — cos e) cos e 

which is satisfied if | < cos0, which is the reason for the choice of = ~. 
Summing up we are left with the condition 

cos 6 — cos e . , 

i < . 13 

cos y — cos e 



We can now finish the proof for (ii) and (iii). If the length y M is reached in step 
2k + 1, then we have 

||w 2 fc+i|| 2 = l + 2fc(cos(0 + e) -cosO + yu)) + 2k 2 (l - cos(> - e)) > M. 

Solving the quadratic inequality in k and using standard trigonometric identities 
we get 



siir(^ + M±£) + M - 1 - sin(0 + ^±£) 
2 sin • 



k > ~ • (14) 



2 

Putting together ( !T3|) and ffl4|) we get 



cos g _ cos e > y sin 2 (0 + a±£) + M - 1 - sin(0 + g±£) 
cos /i — cos e — 2 sin ^ 

Finally we isolate M and arrive at 

(cos e — cos 0) 2 2(cos e — cos 0) sin(0 + ^jjr) 

M < — + — — — — + 1 

sin 2 ^ + sin^ + 

For (ii) replace fx by 3e, for (iii) set e = 0. In both cases the right-hand side 
becomes arbitrarily large when e resp. \x approaches zero. □ 
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